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Abstract 



We study the effects of thermal fluctuations on a small elastic ring. We derive 
analytical expressions for the correlation functions of the Euler angles, for the real 

-6 

space two-point correlation functions and for the probability distribution of writhe, 

o . 

as a function of the persistence lengths that determine the rigidity with respect 
to bending and twist. Fluctuation amplitudes diverge in the limit of vanishing twist 

> 

. rigidity. There is crossover from a small scale regime in which twist and writhe modes 

are uncoupled, to a large scale regime in which twist affects the spatial configurations 
of the ring. 

O . 

Small circular loops of extrachromosomal DNA (plasmids) play an important role in 
biological processes such as gene transfer between bacteria and in biothechnological appli- 
cations where they are used as vectors for DNA cloningtl. The simplest model that captures 
both the topology and the physical properties of such an object is that of an elastic ring, and 
^ . several studies of the writhing instability of twisted rings based on this model, were reported 
in recent yearsSi. However, these studies focused on the elastic response of the ring to me- 
chanical forces and did not consider the effects of thermal fluctuations. These fluctuations 
dominate the physics of macromoiecules and determine all their statistical properties, such 
as radii of gyration, dynamics in solutioni, kinetics of loop formation, and dissociation of 
short DNA segments^ and molecular beaconJI Recently, we developed a theory of fluc- 
tuating elastic filaments, with arbitrary spontaneous curvature, torsion and twist in their 



stress-free state!. Since topological constraints were not taken into account, our analysis was 
limited to linear filaments and could not be directly applied to the study of closed objects 
that have the topology of a ring. In the present letter we present the statistical mechanics 
of weakly fluctuating rings. We calculate the correlation functions of the Euler angles, and 
use them to derive the real space correlation functions, and the complete probability distri- 
bution function of writhe fluctuations. This information allows us to understand the effects 
of twist and bending rigidities on the statistical properties of fluctuating elastic rings. 

Consider an undeformed circular ring with a centerline that forms a circle of radius r in 
the xy plane. The Euler angles that describe this equilibrium configuration are 9q = 7r/2, 
(fo = s and ipo = 0. Here and in the following, the dimensionless contour length s is 
measured in units of r and thus < s < 2tc. In the presence of thermal fluctuations, the 
Euler angles {77} = {9,ip,ip} deviate from their equilibrium values, {770} = {9 , <p , i/j }, and 
the instantaneous deviations are denoted by Sr](s) = r](s) — T)o(s). The topology of the 
ring imposes periodic boundary conditions on the fluctuations of the Euler angles 56(2tt) = 
66(0), 5ip (2ir) = 5ip(0) and <5<^(27r) = 5ip(0), as well as on the fluctuations of the ring 
in three dimensional space, <5x(27r) = <5x(0). Introducing the Fourier transforms, Srj(s) = 
J2nV( n ) eins y where the sum goes over all positive and negative integers and fj(—n) = fj*(n)), 
the above boundary conditions are expressed as 6(0) = <£>(1) = 0. From the general expression 
for the elastic energy of a filament with arbitrary spontaneous curvature, torsion and twist!!, 
one can derive the elastic energy U for small deviations of the Euler angles from their 
equilibrium values, 
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where T is the temperature and a& and a t are dimensionless persistence lengths (given in 
units of the radius of the ring r) associated with bending and twist, respectively, and are pro- 
portional to the elastic moduli that determine the rigidity with respect to the corresponding 
deformations. The quadratic approximation for the energy, Eq. (1), holds in the limit of 



weak fluctuations of the Euler angles, i.e., when their characteristic fluctuation amplitude 
satisfies Srj char <C 1 (in radians). This energy do not depend on modes = —iO(l) and 
<p(0) that correspond to rigid-body rotation of the entire ring, with respect to axes lying in 
the plane of the ring and normal to it, respectively. The expression in the square brackets in 
Eq. (1) contains terms ijpin), 9{n)) that are not multiplied by n and consequently depend on 
the values of the corresponding angles in the stress-free state (the energy depends explicitly 
on the spontaneous curvature of the filament). 

The energy, Eq. (1), is the sum of contributions each of which is a quadratic form that can 
be represented (for n > 2) as a 3 x 3 matrix in the space spanned by fj(n). Since this matrix 
is diagonal in <p(n), with eigenvalues Xi(n) = a^n 2 , the fluctuations of the angle ip decouple 
from those of the angles 9 and ip. The energy of this mode depends on the bending modulus 
only and therefore Sip(s) describes pure bending fluctuations in the xy plane. The remaining 
2x2 matrix in space can be diagonalized by solving for the roots of a quadratic 

characteristic equation, A 2j3 (n) = (a b + a t ) (n 2 + 1) /2=f (a b — a t ) (n 2 + 1) /4 + An 2 a b a t 

We proceed to calculate the correlation functions of the Euler angles, 

(St, (s) Srj' (0)> = £ e ms (rj(n)fj' (-n)) = £ e -*M^, (2) 

n,k K\ k (n) 

where fj(n) and r)'(n) are components of the vector j#(n), <p(n), ^(n)| . The second equality 
in the above equation was derived by expanding the modes fj(n) in the basis of orthonor- 
mal eigenvectors i] k (n) of the quadratic form, Eq. (1), 77(71) = J2k c k(n)r)k(n), and using 
the equipartition theorem to calculate the correlator of the coefficients in this expansion, 
(cfc(n)cfc/(— n')) = [n\k(n)]~ l 5 nn i5kk' ■ The modes with n = 0, ±1 should be considered 
separately, because some of the corresponding eigenvalues vanish. Since the mode Sip de- 
couples from 59 and Sip, the cross-correlators of Sip with 59 and Sip vanish identically, 
(59(s)5ip (0)) = (5ip(s)5ip (0)) = 0. A straightforward calculation yields the following ex- 
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pressions for all the other correlators, 



ir{S<p(s)S<p (0)) 
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(3) 



where the functions fi through / 4 are defined by 



/i(s) = [k(s) + 1/16] coss - l(s) 



f 2 {s) = 2k(s) — cos s, 



(4) 



f 3 ( s ) = [k(s) - 3/16] cos s + l(s) - 1/2, / 4 (s) = [ife(s) + 1/16] sin s. 



with k(s) = (tt — s) 2 /8 — 7r 2 /24, and Z(s) = [(n — s) /4] sins. Inspection of Eqs. (3) shows 



that bare persistence length associated with twist rigidity, a t , plays a fundamentally im- 
portant role: the fluctuations of 59 and of Sip and the cross-correlations between these 
fluctuation modes, diverge in the limit a t — > 0. Therefore, simplified theories of elastic fila- 
ments with non-vanishing spontaneous curvature that do not take into account twist rigidity, 
can not describe the fluctuations of an elastic ring. This problem does not arise in the case 
of straight filaments (with no spontaneous curvature), that can be successfully described by 
a wormlike chain model, in which only bending rigidity is accounted for (a t = 0). 

In Fig. 1 we plot the correlation functions of the Euler angles vs. the dimensionless 
contour distance s, in the interval < s < 2n (for a b = a t = 10). Although the diagonal 
correlations (99, iptp and ipip) decay rapidly with x in the range s < 1, the decay is not 
monotonic and maxima appear at s = 0, n and 2n (these functions are symmetric with 
respect to reflection about s — tt). For 99 and ipip correlation functions, in between these 
maxima there are intervals where fluctuations of Euler angles at points s and s' are in op- 
position, and nodes where these fluctuations are not correlated. This behavior is a direct 
consequence of the ring topology, and the oscillatory patterns can be thought of as standing 
waves produced by interference of two wave packets propagating along two opposite direc- 
tions along the ring. The nondiagonal correlation function 9ip is antisymmetric with respect 
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to reflection about x = it. Surprisingly, Qip correlations vanish in the limit s — > 0, i.e., the 
fluctuations of twist and out-of-plane deviations at two neighboring points are uncorrelated. 
The underlying physics is simple: a short segment of the ring behaves as a straight rod, for 
which twist about the centerline is independent of its orientation with respect to the xy 
plane. 

The real space two-point correlation function, ^[x(s) — x (0)] 2 ^>, can be expressed in 
terms of the correlation functions of the Euler angles. For small s the correlator depends 
only on the bending persistence length, and coincides with that of a straight rod described 
by the wormlike chain model. In Fig. 2 we plot ^[x(s) — x(0)] 2 ^> as a function of s, for 
several choices of the bending and twist rigidities. As expected, the mean square separation 
increases parabolically with s (for small s) and exhibits a maximum at s = irr (determined 
by the geometry of the undeformed ring). Fluctuations suppress this maximum in a way 
that depends on the rigidity parameters: decreasing the twist rigidity a t has a much weaker 
effect than decreasing the bending rigidity a b . 

We proceed to discuss the twist and writhe fluctuations of the ring. The twist (Tw) 
associated with a particular configuration of the ring, can be expressed in terms of the Euler 
angles, 



where 0*3 (s) is the "angular velocity" of twist about the tangent vector^. The writhe (Wr) 
of a given configuration of the ring is defined by0, 



where x and • denote vector and scalar products, respectively, t is a unit tangent to the 
ring contour, and e 2 is a unit vector in z— direction. Expressing t in terms of the Euler 
angles^ in the above expression yields 






dt 



(6) 




(7) 
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Inspection of Eqs. ([5]), (|7|) shows that the total rotation is characterized by the topo- 
logically conserved number, Lk = Tw + Wr that does not depend on the conformation of 
the ring0. In equilibrium, Eqs. (|) and (|7|) give Tw eq = and Wr eq = 1 and although 
fluctuations lead to deviations from these values, the total angle of rotation is conserved, 
2nLk = 2irLk eq = 2ir. Expanding the Euler angles inside the integrals in the deviations 
from their equilibrium values, we obtain the following expressions for writhe and twist fluc- 
tuations, 

SWr = -5Tw = [ 2nr ^S6ds = V in<p{n)6(-n), (8) 

Z7T Jo as „ 

where we used § d5(p = § 56ds = 0. Since both 5ip and 56 enter Eq. SWr vanishes 
both when the fluctuations are confined to the plane of the ring (86 = 0) and when they are 
normal to it (5<p = 0). 

Since the average (if(n)6(— n)) is an even function of n and the sum in Eq. (|8|) goes 
over both positive and negative values of n, we find that (5Wr) = 0, i.e, fluctuations do 
not affect mean writhe. The dispersion of the writhe is given by (5Wr 2 ) = J2n=£o,±i Wr 2 (n), 

1/2 



where Wr(n) = n 



is the rms amplitude of 



6(n)6(-n)) (0(n)fp(-n)) - {(p{n)6{-n) j 
writhe fluctuations at wavelength r/n. Using Eqs. (3) for the correlation functions of the 
Euler angles, we find 



For small wavelengths, n > 1, the amplitude of writhe fluctuations depends only on the 
bending persistence length ab, similarly to the behavior of filaments that fluctuate about a 
straight stress-free configuration@. Indeed, on sufficiently small scales, the ring topology is 
unimportnt and the filament behaves as a straight rod whose ends can freely rotate about 
the tangent direction. Since for rod with zero spontaneous curvature, twist and writhe 
are not coupled through energy (via the dependence of the elastic energy on spontaneous 
curvature, see Eq. (1)) or through topology, writhe fluctuations are independent of twist. 
The crossover to the long wavelength regime at which twist and writhe become coupled, 
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takes place at a length scale £ t = rya t /ab and, therefore, such a regime exists in a ring 
of radius r only if at jab < 1. The straight rod case follows from the above expression by 
substituting r = oo, and since £ t diverges in this limit, we conclude that writhe is not affected 
by twist, independent of the magnitude of the twist rigidity, at. 

Substituting Eq. (P) back into the expression for (5Wr 2 ) yields (6Wr 2 ) = 0.179a 6 ~ 2 + 
0.027a^ 1 a < r 1 - Recall that persistence lengths are scaled by the radius of the ring and, there- 
fore, (5Wr 2 ) ~ r 2 , in agreement with the scaling estimates in referenceEl. Indeed, writhe 



is a quadratic form of Sep and 56, each of which has typical fluctuations of yr/a, where a 
is a characteristic persistent length and, therefore, the characteristic amplitude of writhe 
fluctuations is 5Wr ~r/a. 

We proceed to calculate the probability distribution function of writhe, P (5Wr) = 

5 



5Wr + £ 

n ^o,±i inip(n)6(—n) ) . Using the exponential representation of the ^-function, 



and performing the resulting integrations, a straightforward calculation gives 



OO 

P(^r) = EvrH^- R exp 



\8Wr\ 
Wr(n) 



(10) 



where 



nn 2 (n 2 + a) (1 + a) y/a (n 2 - 2) a b - a t 

7r(n) = (-1) — — - 2 . and a = — . (11) 

2 (n z — 1) sinn (ir^a) a b + n z a t 

The above expression for 7r(n) can be used to calculate all even moments of writhe fluc- 
tuations (odd moments vanish because of the radial symmetry of the undeformed ring), 
8Wr k "j = k\J2% D =2 7T ( n )Wr k (n). The writhe distribution function can be written in the 
form P (nabSWr, ab/dt) and in Fig. 3 we plot this function vs. na^Wr, for ab/a t = 0.1, 5 
and 20. As expected, the probability of large-amplitude writhe fluctuations (large |<5H^r |) 
decreases with increasing twist rigidity (decreasing ab/a t ). The width of the distribution 
decreases sharply when cib/a t is of order unity. With further increase of twist rigidity, the 
narrowing becomes less pronounced and finally saturates for ab/at < 0.1. The distribution is 
non-Gaussian, with an exponential tail at nabSWr ^> 1. The shape of the curve is qualita- 
tively similar to recent computer simulations0, but detailed comparison could not be made 



since the simulations were carried out in the SWr ^> 1 regime (our study is limited to the 
SWr <C 1 range). 

In this letter we presented the statistical mechanics of fluctuating rings. We derived 
analytical expressions for the two-point correlation functions of the Euler angles, for the 
real space correlation functions and for the probability distribution function of the writhe, 
as a function of the persistence lengths associated with the bending and twist deformations of 
the ring. We found that the amplitudes of fluctuations diverge in the limit of vanishing twist 
rigidity and, therefore, theories in which only bending rigidity is taken into account, can not 
be used to model fluctuating rings. We found that a crossover length scale £ t = r^Ja t /ctb 
exists, below which writhe and twist fluctuations are decoupled, and above which twist affects 
the three-dimensional configurations of the centerline of the ring. It appears plausible that 
this crossover does not depend on the topology of the ring and is characteristic of filaments 
with spontaneous curvature in their stress-free state. 

Finally, we would like to comment on the limitations of the present work. Our theory is 
limited to the weak fluctuation regime, in the sense that the deviations of the Euler angles 
from their equilibrium values must be sufficiently small. In order to ensure this, the radius of 
the ring has to be smaller than all the persistence lengths. For simplicity, we only considered 
the case of a ring with a circular cross section. This assumption is not essential and the 
study of a ring with asymmetric cross section and three independent elastic constants, will 
be reported in an expanded version of this letter, together with the study of linear response 
to externally applied moments and forces, based on the fluctuation-dissipation theorem that 
relates the correlation functions to the appropriate response functions!^. 

YR acknowledges support by a grant from the Israel Science Foundation. 



8 



REFERENCES 

1. Electronic address: yr@rabinws.ph.biu.ac.il 

2. Permanent address: Theoretical Department, Lebedev Physics Institute, Russian 
Academy of Science, Moscow 117924, Russia, electronic address: panyukov@lpi.ac.ru 

3. The Encyclopedia of Molecular Biology, Ed. J. Kendrew (Blackwell, Oxford, 1994). 

4. I. Tobias, B.C. Coleman, and M. Lembo, J. Chem. Phys. 105, 2517 (1996). 

5. H. Qian and J.H. White, J. Biomol. Struct, and Dyn. 16, 663 (1998). 

6. P.-G. deGennes, Scaling Methods in Polymer Physics (Cornell Univeristy Press, Ithaca, 
1979). 

7. A. A. Podtelezhnikov and A.V. Vologodskii, Macromolecules 33, 2767 (2000). 

8. N.L. Goddard, G. Bonnet, O. Krichevsky, and A. Libchaber, Phys. Rev. Lett. 85, 2400 
(2000). 

9. S. Panyukov, Y. Rabin Phys. Rev. Lett. 85, 2404 (2000); Phys. Rev. E, 62,. ..(2000). 

10. J.H. White Am. J. Math. 91, 693 (1969). 

11. F.B. Fuller, Proc. Natl. Acad. Set. 68, 815 (1971); Proc. Natl. Acad. Set. 75, 3557 
(1975). 

12. C. Bouchiat and M. Mezard, Phys. Rev. Lett. 80, 1556 (1998). 

13. A.C. Maggs, |cond-mat/9912119| ; |cond-mat/0009182 . 



14. D.A. Beard and T. Schlick, J. Chem. Phys. 112, 7323 (2000). 

15. S. Panyukov and Y. Rabin, in preparation. 



9 



Figure captions 

Fig. 1: Plots of two-point correlation functions of Euler angles vs. the dimensionless 
contour distance between the points, s, in the interval < s < 2ir: (59(s)59(0)) (box), 
{5(p(s)8(f(0)) (diamond), (5i/j(s)8iIj(0)) (cross) and {5 9 (s) 5^(0)) (solid line). The bare per- 
sistence lengths are a b = a t = 10. 

Fig. 2: Plot of real space two-point correlation function ^[x(s) — x(0)] 2 ^ vs. the 
dimensionless contour distance between the points, s, in the interval < s < 2n. The 
parameters are = 10, a t = 10 (solid line), a& = 10, a t = 1 (box) and — 1, a t — 10 
(cross). 

Fig. 3: Plot of probability distribution function of writhe vs. ira b 5Wr for a b /a t = 0.1 
(solid line), 5 (cross) and 20 (box). 
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